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Abstract 

NA44 has measured mid-rapidity deuteron spectra 
from A A collisions at ^s„„ « 18 GeV at the CERN 
SPS. Combining these spectra with published p, p 
and d data allows us to calculate, within a coales- 
cence framework, p and p source sizes and phase space 
densities. These results are compared to n"^ source 
sizes measured by Hanbury-Brown Twiss, HBT, in- 
terferometry and phase densites produced by com- 
bining pion spectra and HBT results. We also com- 
pare to pA results and to lower energy (AGS) data. 
The p source is larger than the proton source at 
y'Snn = 17.3 GcV. The phase space densities of 7r+ 
and p are not constant but grow with system size. 
Both TT"*" and proton radii decrease with and in- 
crease with ^Jsnn- Pions and protons do not freeze- 
out independently. The nature of their interaction 
changes as ^/Snn and the n/p ratio increases. 

Relativistic heavy ion collisions provide a mech- 
anism to heat and compress nuclear matter to tem- 
peratures and energy densities comparable to those of 
the early universe when it was still a plasma of quarks 
and gluons. Such a state may be fleetingly restored 
in these collisions where temperatures of T = 168 ± 3 
MeV and energy densities e = 3 GeV/fm'^ are ob- 
served ||l|, ^. These values are close to those of the 
phase transition found in lattice calculations If 
such a hot and dense state were formed one would 
expect a large increase in entropy and possibly a sat- 
uration of the density of particles in phase space. The 
coalescence of nucleons into deuterons is sensitive to 
both their spatial and momentum correlations. In 



this paper we report p, p and 7r+ source sizes mea- 
sured by coalescence and interferometry, and com- 
bine these with single particle spectra to derive phase 
space densities. The phase space densities depend on 
temperature, chemical potentials, and velocity fields 
in the system. This description of the final hadronic 
state serves as a boundary condition for models of 
possible quark gluon plasma production. We vary 
the total size of the system by studying PbPb, SPb, 
SS and pPb collisions. We also compare our results 
to lower energy AGS data where the 7r/p ratio is much 
lower. This comparison shows that the freeze-out of 
pions and protons is coupled. 

NA44 is a focusing spectrometer that uses 
conventional dipole magnets and superconducting 
quadrupoles to analyze the momentum of the pro- 
duced particles and create a magnified image of the 
target in the spectrometer [|[ |, g @, |, §, The 
systematic errors on the deuteron yields range from 
14% for SPb to 9% for PbPb. The p and p spectra 
are corrected for feed-down from A and E decays us- 
ing a GEANT simulation with the (A/p) and (S/p) 
ratios taken from the RQMD model ||, |ll|. The 
systematic error was estimated by varying these ra- 
tios by ±25%. These errors are slightly correlated 
for p and p. Fig. |^ shows NA44 deuteron spectra and 
previously published proton spectra at y = 1.9-2.3 as 
a function of my /A |l^. The centrality is ~ 10% 
for SS, SPb and PbPb. The spectra get flatter for 
the larger systems consistent with a higher temper- 
ture and/or stronger sidewards flow. As expected 
from coalescence, the slopes are similar for protons 
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Figure 1: Invariant cross-sections vs. rriT/A at 
y=1.9-2.3, for protons (a) and deuterons (b) from 
central SS, SPb and PbPb collisions §, |l| 



and deuterons. This was also found for lower energy 
data, 



^ The deuteron inverse slopes (in ttit) and 

yields are listed in Table A comprehensive analysis 
of all NA44's proton and light clusters spectra will be 
given in a later paper. 

The model of deuteron production by final state 
coalescence of protons and neutrons with small rela- 
tive momenta states that the production of deuterons 
with a certain velocity is proportional to the num- 
ber of protons and neutrons that have similar veloc- 
ities |l6|. This model successfully describes mea- 
sured deuteron distributions in intermediate energy 
heavy ion collisions and high energy pA collisions, 
[ p7[ . Near mid-rapidity, direct production of dd pairs 
is small due to the high dd mass threshold of 3.75 
GeV/c^, and pre-existing deuterons are unlikely to 
survive the many collisions required to shift them to 
mid-rapidity. Since coalescence depends on the dis- 
tribution of nucleons, one can determine a nucleon 



System 


Fit Range 

of niT — m 

MeV/c^ 


Inverse 
Slope 

MeV/c^ 


dN/dy 

(10-2) 


SS 
SPb 
PbPb 


0-520 
0-520 
160-520 


320 ± 149 
300 ± 91 
379 ± 13 


39 ± 13 
153 ± 23 
390 ± 20 



Table 1: Deuteron inverse slopes and yields. System- 
atic and statistical errors have been added in quadra- 
ture. The errors are dominated by statistics and the 
extrapolation out of the acceptance. The PbPb fit is 
from H. 



source size from the ratio 

B2{P) = 



dP=> 



(1) 



where the deuteron momentum P is twice the proton 
momentum p p8[ | . Since we do not measure neutrons 
we assume that the spectra have the same shape and 
take the n/p ratio to be 1.06 ± .04 from RQMD 0. 
At ^s„„ = 4.9 GeV the measured n/p ratio is 1.19 ± 
.08 independent of mj- |l^. RQMD is in reasonable 
agreement with this result. 

To facilitate comparison with NA44's pion inter- 
ferometry results, we assume a Gaussian distribution 
of the proton source. If one also assumes a Gaus- 
sian wave-function one can solve for the source size 
analytically PQ] 
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where is the proton mass and (5 = 2.8 fm accounts 
for the size of the deuteron. In reality the deuteron 
wave function is not Gaussian but is more accurately 
represented by the Hulthen form 
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with a = 0.23 fm-i and f3 = 1.61 fm'^ |2^. The 
convolution of such a wave function with a gaussian 
source cannot be done analytically but is straightfor- 
ward numerically to solve for the source radius Rh, 
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Figure 2: Comparison of proton source radii versus 
the coalesence parameter shown by •, is 

derived numerically using the Hulthen wave function 
in Eqn.||. The solid line is a fit to the numerical 
results. The dotted and dashed curves show Rq from 
Eqn|2| with and without a correction for the deuteron 
size. 



H. Note our Rh is the i?o of (H. Figure | shows 
a comparison of Rh and Rq versus the coalesence 
parameter i?2. 

Since Ru is sensitive to both the transverse and 
longitudinal size of the source, when comparing to 
HBT results it is best to compare to (i?^ ' ^|| ) ^ (Eqn. 
6.3 of ^^), where i?x and parametrize the extent 
of the source perpendicular and parallel to the beam 
||||. NA44 has published HBT results in the Pratt- 
Bertsch frame in which the sum of the longitudinal 
momentum of each pion pair is zero. In this scheme 
the radius in the sidewards direction R^ — R± and 
the longitudinal radius Ri ~ i?|| In this paper 

we will therefore compare Rh to ■ Ri)3. These 
parameters can be thought of as "lengths of homo- 
geneity" of the source |2^, One can think of the 
radii as lengths scales of the velocity and/or tem- 



perature gradients. They represent snapshots of the 
hadronic system at freeze-out which may occur at dif- 
ferent times for 7r+ and protons. However since the 
cross sections for tttt, pn and pp collisions are compa- 
rable the freeze-out times should be close. 
A particle's phase space density is defined as 



dp-^dx^ 



(4) 



For a system in chemical equilibrium at a tempera- 
ture T and chemical potential /i 



f{E) = 



1 



,{E-t,)/T ^ I 



(5) 



where E is the energy and ±1 selects bosons or 
fcrmions. For a dilute system, i.e. / <C 1, Eqn. |^ 
gives 

Since Ed = En + Ep, Eqn. ^ implies that 



/d(P,x) = /p(p,x)/„(p,x). 



p 



/p(P,x)2 (7) 



A more general form of this relation was derived in 
Eqn. 3 of assuming only that the system is hot 
and large compared to the deuteron. Averaging fp 
over X gives 
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where the factor of 3 accounts for the spin of the par- 
ticles. For pions NA44 has measured the source size 
in 3 dimensions with HBT, as well as single particle 
spectra. Some of the pions come from the decay of 
long-lived resonances such as r],r]' and uj. These pi- 
ons reduce the strength of the correlation function A, 
which typically is < 1. The fraction of pions which 
do not come from resonances is \/A, This has 

been shown experimentally for e'*~e~ collisions and 
for RQMD simulations of PbPb collisions |3^, |l). 
Dividing ^/Xd^N^^/dp^ by the Lorentz invariant vol- 
ume, [pTl |33| gives 
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Again the factor of 3 accounts for the pion's spin de- 
generacy. Here Ri is the extent of the source along 
the beam direction; Ro the extent in the outward di- 
rection, ie towards the observer and Rg measures the 
source in the sidewards direction, perpendicular to 
both the beam axis and the line of sight to the ob- 
server. The Roi term is the "out-longitudinal" cross 
term [|5|. For PbPb collisions setting Roi = in the 
HBT fit increases (U) by 9% ± 10%, §. For pPb, 
SS and SPb we assume Roi = but add a systematic 
error of 13% to (/^). For pPb deuteron spectra are 
not available and (/p) was calculated using Eqn. ^, 
replacing the last term with — . Rinv was deter- 
mined from pp HBT data [Q. For PbPb Rh and 
Rinv agree within their errors of ~ 5%, [ p4| . 

We can test the usefulness of these coalesence 
methods using RQMD, coupled with a coalesence 
afterburner, |^5j. Figure ^ shows a comparison of 
(/p) and Rh for both data and RQMD for SS and 
PbPb. For the data both (/p) and Rh are larger in 
PbPb but for RQMD only Rh increases from SS and 
PbPb while (/p) stays constant. This invariance of 
(/p) may be an artifact of the coalesence mechanism 
used, which ignores the requirement of a third parti- 
cle. Since (/p) is constant in the model the increase 
in proton multiplicity from SS to PbPb is accommo- 
dated by a large increase in Rh- However, Rh is less 
than the average transverse radius of freezeout, 5.4 
fm for SS and 10.3 fm for PbPb, indicating that coa- 
lesence is not sensitve to the full size of the source. A 
similar situation occurs in pion interferometry where 
correlations between position and momentum cause 
the observer to only "see" the side of the source clos- 
est to her, 1^. Since these correlations get stronger 
as the particles get faster the size of the source drops 
with mr, |29| ||. 

Figure | shows the system dependence of the phase 
space densities and source radii for tt^ , p and p. The 
p and p radii for PbPb are consistent with coalescence 
data at pt = 0, and pp interferometry results |^, Q . 
The TT"*" and p phase space densities generally increase 
with system size. We find that 
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Figure 3: Proton phase space densities (/p) (a),(b) 
and radii Rh (c),(d) versus rtiT for ^s„„ — 17.3 GeV 
for data • and RQMD □. The shaded bands indicate 
the estimated systematic error on the correction for 
weak decays. 



(/p)«(/p) «(/.+ ) <(/.-)«! 
For SPb and PbPb (/7r+) was calculated in |3^ using 
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a similar equation to Eqn. ^. However a parametriza- 
tion of the pion spectrum was used rather than the 
spectrum itself. The authors of concluded that 
the pion phase density was universal at freeze-out 
but this is not the case, since (/^+) is considerably 
smaller for pPb and SS than for SPb and PbPb. 

For pions (i?^ ■ Ri)^ increases steadily with the 
number of participants and for PbPb there is a rapid 
increase in the radii parameters with multiplicity, § . 
At low pT Rh does not change much from SS to 
PbPb (nor with centrality for PbPb) despite the in- 
crease in the proton multiphcity by « 3, see Fig. |l|(a). 
However the dependence of Rh is weaker for 
PbPb than for SS, see Fig ||. The extra protons 
mainly increase the proton phase space density (fp). 

Because of their large annihilation cross-section, 
one might expect that antiprotons (particularly those 
at low pt) would be emitted only from the surface of 
the system and would have a larger RMS freeze-out 
radius than protons. Our data are consistent with 
this idea. An alternative view assumes that pro- 
ton and antiprotons are produced in the same vol- 
ume but that antiprotons are suppressed in the in- 
terior of the source |40 . Applying this idea to our 



data would imply that antiprotons are emitted only 
from within 1.0 ± 0.2 fm of the surface |Q. Recent 
AuPt and AuPb results from E864 at y/s^ = 4.9 
GeV and low pr imply that i?^ = (4.0 ± 0.2) fm, 
R^r = (2.2±0.9±0.6)fmand (fp) = (4.0± 1.9) • lO^^ 
|{41| . At Y^s^ = 4.9 GeV antiprotons are mainly 
produced in primary nucleon-nucleon collisions and 
so they may have a smaller source size than protons. 
Since Rh is a gaussian radius it is necessary to mul- 
tiply it by in order to compare it to a hard sphere 
with the same RMS radius, If this is done the 

antiproton source is roughly equal to the size of the 
colliding Au nuclei at ^/Snn = 4.9 GeV. 

In order to study the energy dependence of freeze- 
out we compare our PbPb data at y/Snn = 17.3 GeV 
to AGS AuAu data at ^^Snn — 4.9 GeV. Since our 
data are not at mid-rapidity but at y=2 wc have com- 
pared results at the same scaled rapidity y — ■^ybeam- 
Figure || shows the phase space densities and source 
radii for PbPb and AuAu collisions as a function of 
rriT- At a given ttit, (/7r+) increases with y/Snn while 
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Figure 4: Phase space densities (/) and source radii 
for 7r+ and p at (pt) ~ 240 MeV/ c, and for p at 
(PT> ~ 490 MeV/c. For pPb, = 30 GeV and 

the proton points are derived from pp HBT data M . 



(fp) decreases. Fitting {fjr+) to Eqn. js] gives ^+ = 0, 
within errors, for both energies while fJ.p/T decreases 
with y/Snn- Since E — mTCOsh{y), Eqn. || also im- 
plies that / be exponential in ttit for / ^ 1. How- 
ever if the system is boosted due to transverse flow, 
f(mT) will become flatter |Q. This effect is propor- 
tional to mass. The data support this scenario since 
the (fp) distributions are much flatter than the (/7r+) 
ones. The mx distribution of (fp) becomes flatter as 
^s„„ increases because of an increase in flow and/or 
freeze-out temperature. However in general the ve- 
locity profile cannot be determined without knowing 
the density profile and so a determination of a mean 
velocity from (fp) is beyond the scope of this work. 
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Several hydrodynamical models have interpreted 
the HBT source radii as "lengths of homogeneity" 
which should decrease with increasing mx and this 
is consistent with our data |2^. Both pion and 
proton radii increase with ^s„„. However (/7r+) in- 
creases with y/Snn while (fp) drops. Since p/p <C 1 at 
both SPS and AGS energies §, ^ 0, we know that 
most protons observed near mid-rapidity are rem- 
nants of the target or projectile that were slowed 
down by multiple collisions. At the higher energy 
the protons occupy a somewhat larger volume and 
they are spread over a larger momentum (i.e. y, ttit) 
range. Therefore {fp) drops with ^Jsnn- 

For pions the situation is different. At ^Jsnn — 4.9 
GeV, they are outnumbered by protons and so their 
freeze-out is driven by that of the nucleons. At y^Snn 
= 17.3 GeV, they are the most numerous particle and 
control freeze-out. Since fjjrTr < o'prr {f-n+) increases 
with y^Snn- Note that the ratio {f-,z+)/{fp) increases 
by a factor of about 16 from ^Snn ~ 4.9 GeV to 17.3 
GeV while the ir/p ratio only increases by a factor of 
7. 

Using the (anti)dcuteron as a measure of the 
nucleon-nucleon correlations we have used a coales- 
cence formalism to make the first measurements of 
p and p source radii and phase space densities as a 
function of ttit at ^s„„ = 17.3 GeV. At y/Snn = 4.9 
GeV the antiproton source is smaller than the pro- 
ton source while at ^/Snn = 17.3 GeV it appears to 
be larger. We have compared the proton data to our 
TT"*" radii and phase space densities derived from HBT 
and single particle results as a function of system size 
and y/Snn- This comparison reveals a linkage between 
proton and pion freeze-out that changes as the 7r/p 
ratio increases. At y/Snn = 4.9 GeV the hadronic 
system is held together by protons while at y/Snn — 
17.3 GeV it is held together by pions. 

We thank the NA44 collaboration for permission 
to publish deuteron spectra and much valuable input. 
Thanks also to G. Bertsch, U. Heinz, S. Mrowczyhski 
and S. Pratt for helpful discussions. 
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Figure 5: Phase space densities (a) and radii (b) for 
7r+ and p versus tot for ^^Snn = 17.3 and 4.9 GeV. 
The shaded bands indicate the estimated systematic 
error on the correction for weak decays. The E866 
points used data from O pi pf. 
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